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ABSTRACT 

We have investigated the variation of magnetic hehcity over a span of several 
days around the times of 11 X-class flares which occurred in seven active regions 
(NOAA 9672, 10030, 10314, 10486, 10564, 10696, and 10720) using the mag- 
netograms taken by the Michelson Doppler Imager (MDI) on board the Solar 
and Heliospheric Observatory {SO HO). As a major result we found that each 
of these major flares was preceded by a significant helicity accumulation, (1.8- 
16) xlO^^ Mx^ over a long period (0.5 to a few days). Another finding is that the 
helicity accumulates at a nearly constant rate, (4.5-48) xlO^'' Mx^ hr~^, and then 
becomes nearly constant before the flares. This led us to distinguish the helicity 
variation into two phases: a phase of monotonically increasing helicity and the 
following phase of relatively constant helicity. As expected, the amount of helic- 
ity accumulated shows a modest correlation with time-integrated soft X-ray flux 
during flares. However, the average helicity change rate in the first phase shows 
even stronger correlation with the time-integrated soft X-ray flux. We discuss 
the physical implications of this result and the possibility that this characteristic 
helicity variation pattern can be used as an early warning sign for solar eruptions. 
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1. INTRODUCTION 

Magnetic helicity is a measure of twists, kinks, and inter-linkages of magnetic field lines 
(Berger & Field 1984) and has been an important parameter in solar dynamo theories (Parker 
1955). While the somxe of the magnetic helicity lies below the surface of the Sun, it was 
recently recognized as a useful parameter in describing solar features observed above surface 
such as spiral patterns of sunspot fibrils, helical patterns in filaments and coronal mass 
ejections (CMEs; for a review, see Brown et al. 1999). Naturally magnetic hehcity studies 
have been directed to the energy buildup and instability leading to eruptions and CMEs (e.g.. 
Rust 2001; Kusano et al. 2004; Phillips et al. 2005). More recently, several studies were 
carried out to relate the change of magnetic helicity to the problem of impending or triggering 
solar flares. Moon et al. (2002a) studied the magnetic helicity change around major flares 
to flnd its rapid helicity change before flares, and concluded that a sudden helicity injection 
may trigger flares. Moon et al. (2002b) applied the same approach to seven homologous 
fiares in the active region, NCAA 8100, over a period of 6.5 hr to find a good correlation 
between the amount of incremental helicity and the soft X-ray fiux during each homologous 
fiare. The results from both studies thus point to the idea that the hehcity change occurring 
over short timescales (around a half hour) can be a significant factor in triggering fiares. 
Kusano et al. (2003) proposed annihilation of magnetic helicity as a triggering mechanism 
for solar flares. Numerical simulations were carried, which show that, if the helicity is sharply 
reversed within a magnetic arcade, reconnection quickly grows in the helicity inversion layer, 
driving explosive dynamics. Yokoyama et al. (2003) studied flare activities in active region 
NCAA 8100 to find that most of the fiare events occurred about half a day after the helicity 
injection rate changed its sign, and the positions of Ha emission in fiares well correspond 
to the hehcity inversion fines in space. Sakurai & Hagino (2003) studied two active regions 
appeared in 2001 (NCAA 9415 and 9661), both of which have produced X-class fiares. Their 
finding was, on the contrary, that the magnetic helicity integrated over the regions evolved 
slowly and did not show abrupt changes at the time of the flares, although the distributions 
of magnetic helicity changed significantly over a few days in the regions. 

In this paper, we study long term (a few days) variations of the magnetic helicity around 
major X-class flares. While some of the above studies suggest short-term helicity change as 
an important topic for flare triggering, Hartkorn & Wang (2004) found that the rapid helicity 
change at the time of a flare can occur as an artifact under the influence of flare emission 
on the spectral line adopted in MDI measurements. This means that a short term variation 
during strong fiares can hardly be measured with enough accuracy. However we can study a 
long-term variation of the magnetic helicity with active region magnetograms gathered over 
many days excluding the times of fiares. Our research is therefore focused on a possible 
characteristic helicity evolution pattern that is associated with fiare impending mechanisms. 
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2. 



CALCULATION OF MAGNETIC HELICITY 



By magnetic helicity, we refer to the relative magnetic helicity in the rest of this pa- 
per, i.e., the hehcity relative to that of the potential field state. With the time-dependent 
measurement of longitudinal magnetic fields in the photosphere, we can only approximately 
determine the change rate of the relative magnetic helicity (Demouhn & Berger 2003). We 
further use a simplified expression for the helicity change rate (Chae 2001) given by 



where is the normal component of the magnetic field; Ap is the vector potential of 
the potential field; vlct represents the apparent horizontal motion of field lines; dS is the 
surface integral element and the integration is over the entire area of the target active region. 
Although this expression does not explicitly include the helicity change by the vertical motion 
of field lines (see Kusano et al. 2002), Demoulin & Berger (2003) pointed out that it actually 
accommodates both the vertical and horizontal motions of fiux tubes as far as no fiux tube 
newly emerges from or totally submerges into the surface. As we will show later in presenting 
the result of helicity calculation, this requirement is not always met. 

We determine the quantities in the above equation using full disk MDI (Scherrer et al. 
1995) magnetograms following the procedure described in Chae & Jeong (2005). First, we ap- 
proximately determine Bn from the line-of-sight magnetic field Bi in the MDI magnetograms, 
simply considering the projection effect, i.e., Bi — BnCOSip where ip is the heliocentric angle 
of the point of interest, assuming that the magnetic field on the solar photosphere is normal 
to the solar surface. Second, Ap is calculated from Bn by using the Fast Fourier Trans- 
form method as usual. The extent of the spatial domain of the Fourier transform is taken 
about twice wider than the active region area in order to minimize the artifacts arising from 
the periodic boundary condition in the fast Fourier transform (Alissandrakis 1981). Third, 
vlct is calculated using the local correlation tracking (LCT) technique (November & Simon 
1988). For local correlation tracking, we ahgn all magnetograms in each event to the first 
image of the data set after correcting the differential rotation. We set the FWHM of the 
apodizing window function to 10" and the time interval between two frames to 60 minutes, 
and performed LCT for all pixels with an absolute flux density greater than 5 G. Only the 
pixels with cross correlation above 0.9 are considered. 

In selecting data we found that use of 1 minute cadence full-disk MDI (Scherrer et 
al. 1995) magnetograms is adequate for our purpose of investigating the long-term helicity 
evolution. However, there are occasionally found data gaps in the 60 minute cadence data 
set in which case we supplement the data gaps with 96 minute MDI magnetograms. The 
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time interval of the supplemented data set is therefore not longer than 96 minute. To reduce 
the effect of the geometrical projection, we selected the active regions lying within 60% 
of the solar radius from the apparent disk center. Note that we use only full disk MDI 
magnctograms that have 2"x2" pixel size. Therefore the LCT velocities calculated here may 
have been systematically underestimated compared with the LCT velocities calculated with 
the higher resolution (0.6"x0.6") MDI data (Longcope et al. 2007). 

After the helicity change rate is determined as a function of time, we integrate it with 
respect to time to determine the amount of helicity accumulation: 



where to and t arc the start and end time of the helicity accumulation, respectively. If to is 
a time when the magnetic field is in the potential state, AH is simply the helicity, H{t), at 
time t. However, there is no guarantee that we can observe, by chance, an active region in 
the potential energy state. We therefore set to as the earliest time without significant helicity 
accumulation at which the average value of the helicity change rate over four hours is less 
than our nominal threshold in helicity change rate, 1x10^° Mx^ hr~^. If that time cannot 
be determined, we define to as the time when the data set starts or when the previously 
accumulated helicity is released by a flare. The exact time of to here is unimportant because 
it is only a trial value. After determining H{t), we redeflne to as the time when the resulting 
helicity starts to increase from a nearly constant value. 



Wc present the helicity variation calculated for seven active regions in Figures 1 and 
2. In both figures we plot the magnetic helicity accumulation together with the GOES soft 
X-ray light curve and magnetic fiux as functions of time. The soft X-ray light curve is shown 
to indicate the fiare times and the magnetic fiux is shown to check the above-mentioned 
requirement for the approximation made in Equation (1). Note that the fluxes shown in this 
paper are total unsigned magnetic flux, i.e., sum of the absolute amounts of positive and 
negative fluxes, because net magnetic flux may show little change despite signiflcant flux 
change in each polarity. 

For the events shown in Figure 1 we can see that the helicity accumulates at a monotonic 
rate of change about 0.5-2 days before the flare onset, and then becomes almost constant 
before the flares. For convenience, we distinguish the magnetic helicity variation in two 
stages: a phase of monotonically increasing helicity (phase I) and the following phase of 
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relatively constant helicity (phase II). This pattern is obvious for the four flares (2001 October 
25, 2004 November 7, and 2005 January 16 and 17). For the 2005 January 15 event, the 
helicity increased up to 22:00 UT on January 14 and then decreased afterward. In this case, 
we do not consider that the flare occurred in phase II. It is then noted that these flares took 
place after a signiflcant amount, ^(1.8-11) xlO^^ Mx^, of helicity accumulation. 

In Figure 2, we show the result for the other four active regions. Like the events in 
Figure 1, these flares also occurred after a signiflcant helicity accumulation, ~(1.9-16) xlO"^^ 
Mx^. However, they occurred in the middle of the continuous helicity accumulation, unlike 
those events shown in Figure 1. In other word, the flares in Figure 2 occurred in phase 
I, while those in Figure 1 occurred in phase II. One common trend is, however, that all 
the events are apparently associated with a considerable amount of helicity build-up before 
the flares, whether they occurred in phase I or in phase II. In case where flares occur in 
phase II, it may imply that solar active regions can wait for major flares after the helicity 
accumulated to some limiting amount. This is seemingly contrary to the general belief that 
a flare occurs as soon as the system reaches some threshold. An active region may evolve to 
a certain stage where the helicity no longer increases, and the system waits until it unleashes 
the stored energy by producing flares due to certain mechanism of triggering. 

Since we claimed that these large flares are always preceded by signiflcant accumulation 
of helicity, as a reference we check the corresponding helicity variation in non-flaring times. 
We show such data in Figure 3. For all active regions under investigation, the amount of 
helicity change during nonflaring periods (Fig. 3) is much less than that around the major 
flare time. This convinces us that the above monotonically increasing helicity before major 
flares is a process associated with the flares and is not occurring in nonflare times. Another 
point to note in Figure 3 is that not only the helicity but also the total unsigned magnetic 
flux changes much less during the nonflaring time compared with the period before major 
flares. This implies that the increase of the magnetic helicity before major flares is, in part, 
related to the simultaneous increase of total unsigned magnetic flux. 

It is also worthwhile to mention how the characteristic pattern of the helicity variation 
found here will depend on the sign of helicity. In our result obtained for seven active regions, 
similar amounts of both of positive and negative helicity were accumulated continuously and 
simultaneously during the whole time. It is therefore unlikely that counting the helicity in 
one and the other polarity separately yields a signiflcantly different conclusion. On the other 
hand, some studies suggested that the sign-reversal of the helicity injection rate is important 
for flare activity and we need to compared them with the present result. Kusano et al. (2003) 
emphasized spatially sharp reversal of helicity sign triggers magnetic reconnection based on 
model simulation. While the model prediction is interesting and compelling, we excluded 
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from the outset (see §1) the rapid hehcity change during the flare time due to observational 
hmitation. Our conclusion is valid only for the long term variation of helicity. Yokoyama 
et al. (2003) have found that flares tend to occur after reversal of helicity injection rate 
changed its sign. Although this is occasionally seen in our samples, (i.e., in the case of AR 
10030 and AR 10720) as well, it is not always the case and we are unsure whether this is a 
necessary condition for the flares. More often than not, the helicity either remains constant 
or increases in one sign when the flare occurs. 

4. Correlation with Soft X-ray Flux 

Wc compare the helicity change rate, dH/dt and accumulation amount, AH with the 
time integrated soft X-ray flux taken as the proxy for the flare energy release. In addition wc 
check the helicity accumulation time, r, deflned as the time interval of helicity accumulation 
measured from and the first coming flare. Prior to make such a comparison, the range of 
uncertainty of each quantity needs to be known. In general it is hard to trace all the possible 
uncertainties involved with each quantities in Equation (1). Fortunately, our targeted quan- 
tity given by Equation (2) involves integration in space and time and the uncertainty in each 
measured quantity is not propagating, but rather may cancel out in the process of spatial 
and time integration if it is random in nature. We thus focus the uncertainty estimate only 
on the linear approximation of the helicity variation that we arc after. Wc first find out the 
best-fit linear function to the points AH{ti) lying in phase I (i.e., to <ti < to + r) in the form 
of F{t) = a{t — io) + F{tQ). Next we calculate the standard deviation, a, of the scatter points 
with respect to this linear function, and plot two additional lines corresponding to the ±0" 
levels of the scatter points. Finally, we read the y-axis offsets of these two fines to determine 
cah and the x-axis offsets to determine cr^, respectively. In addition, we calculate the uncer- 
tainty of the slope a itself in the form of {AH — a ah) + cTt) < a < {AH + a ah) /{t — Ct)- 
The center value of a here is taken as the average helicity change rate in the rest of this 
paper. Therefore, by the average helicity change rate, we do not mean the average of the 
quantities given in Equation (1), but wc refer to the best fit slope to the helicity variation 
(eq. [2]) in phase I. The uncertainties shown in Figure 4 and Table 1 are those associated 
with our linear function fit only. 

In Figure 4, we plot, as symbols, the helicity parameters against the GOES soft X- 
ray fiuxcs integrated over the fiaring time {Fx, hereafter). Each symbol is identified with 
the event ID number in the figure together with uncertainty range represented by the bar 
(see also Table I). The solid fines show the least-squares linear fits to the data points. The 
correlation coefficients (CCs) of the finear fits are also given in each panel. Figure 4a shows 
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that there is a fairly good correlation (CC=0.86) between the hehcity change rate and Fx- 
The amount of helicity accumulation also shows a modest correlation with Fx (CC=0.68) 
as shown in Figure 4b, although not as good as for the helicity change rate. On the other 
hand, the correlation between helicity accumulation time r and the soft X-ray flux is very 
poor with a weak tendency that the longer accumulation time r, the weaker soft X-ray flux 
Fx (Fig. 4c). 

We initially expected, on a general basis, that the helicity change AH would strongly 
correlate with Fx- It is therefore puzzling why the helicity change rate dH/dt shows even 
a better correlation with Fx in Figure 4. As a possibility, wc considered that r may be a 
factor in complicating the relationship between Ai/ and Fx- A intriguing idea is that the 
magnetic energy decays much faster than the magnetic helicity in the presence of magnetic 
diffusion (Berger 1999). We thus compare AH with r in Figure 4d, which unfortunately 
shows no obvious correlation between them. The small number of events used in this study 
is another restriction for finding a trend here. With the present result alone, it is fair to 
presume that the weaker correlation between Fx and AH may arise from our inaccurate 
determination of the helicity accumulation amount due to unknown initial time of helicity 
build-up. 



5. SUMMARY 

We have investigated the variation of magnetic helicity over a time span of several days 
around the times of 11 X-class flares which occurred in seven active regions using MDI 
magnetograms to find the following results. 

First, a substantial amount of helicity accumulation is found before the flare in all the 
events. The helicity increases at a nearly constant rate, (4.5-48) x 10^° Mx^ hr~^, over a 
period of 0.6 to a few days, resulting in total amount of helicity accumulation in the range 
of (1.8-16) X 10'^^ Mx^. Such a wide range of helicity accumulation indicates that each active 
region has its own limit of helicity storage to keep a stable magnetic structure in the corona. 
The finding of a monotonically increasing phase is similar to the earlier one by Sakurai & 
Hagino (2003) that the magnetic helicity integrated over the regions evolved slowly and did 
not show abrupt changes at the time of the flares. The helicity increase over days before 
the flares reconfirms the conventional idea that helicity accumulation by a certain amount 
is necessary for a large fiare to occur (Kusano et al. 1995; Choe & Lee 1996). 

Second, there is a strong positive correlation between the average helicity change rate 
of phase I and the corresponding GOES X-ray fiux integrated over the fiaring time. The 



-8- 



amount of helicity accumulation during phase I also correlates with the soft X-ray flux, as 
expected, but the correlation is stronger with the helicity change rate. This result probably 
implies that the helicity change rate is more accurately determined than the amount of 
helicity change itself as the initial time of helicity build-up is poorly determined. 

If the above correlations hold for a large number of events, we may predict the flare 
strength (the integrated X-ray flux) based on the helicity change rate. Monitoring of helicity 
variation in target active regions may also aid the forecasting of flares. A warning sign of 
flares can be given by the presence of a phase of monotonically increasing helicity, as we 
found that all the major flares occur after signiflcant helicity accumulation. As a reference 
we have checked helicity variation of the six active regions in non-flaring times to find much 
lower helicity change rates compared with those around the major flares. We thus conclude 
that the relative magnetic helicity can be a powerful tool for predicting major flares. 
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Tabic 1. List of flares, hclicity and acciimiilation time 



ID 




Flares 


AR 


Peak time 




\dH/(U}' 1 












number 


(C/T) 


(10-1 J/m2) 


(10«Ma:2/ir--i) 


(lO'i^Mx^) 


ihr) 


1 


X 


1.3 on Oct 25, 2001 


9672 


15:02 


2.3 


6.2±0.9 


1.8±0.1 


29±2 


2 


X 


3.0 on Jul 15, 2002 


0030 


20:08 


1.4 


13.3±1.5 


1.9±0.1 


14±1 


3 


X 


1.5 on Mar 17, 2003 


0314 


19:05 


1.3 


6.4±1.6 


3.9±0.5 


63±8 


4 


X 


1.5 on Mar 18, 2003 


0314 


12:08 


1.3 


13.3±1.4 


2.3±0.1 


17±1 


5 


X 


18 on Oct 28, 2003 


0486 


11:10 


20.0 


48.4±6.3 


14.5±0.9 


30±2 


6 


X 


10 on Oct 29, 2003 


0486 


20:49 


9.1 


46.8±3.0 


15.9±0.5 


34±1 


7 


X 


1.2 on Feb 26, 2004 


0564 


02:03 


0.75 


4.5±0.4 


3.1±0.1 


70±3 


8 


X 


2.2 on Nov 07, 2004 


0696 


16:06 


2.1 


19.8±0.7 


10.7±0.2 


54±1 


9 


X 


1.3 on ,Jan 15, 2005 


0720 


00:43 


1.3 


22.2±4.1 


4.2±0.4 


19±2 


10 


X 


2.8 on .Jan 15, 2005 


0720 


23:00 


6.6 


22.5±1.2 


4.3±0.1 


19±1 


11 


X 


4.1 on Jan 17, 2005 


0720 


09:52 


9.1 


40.8±5.5 


3.2±0.2 


8±1 



^-Integrated GOES X-ray flux. 

''Average helicity change rate of phase I. 

'^The amount of heUcity accumulation during phase I. 

"iHelicity accumulation time. 
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Fig. 1. — Time variations of helicity accumulation, magnetic flux, and GOES X-ray flux for 
three active regions. The helicity is shown as cross symbols and the magnetic flux is shown 
as diamonds. The GOES X-ray flux is shown as the dotted lines. Phase I, the interval over 
which the helicity accumulation is considered, and phase II, the following phase of relatively 
constant helicity, are marked. 
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Fig. 3. — Time variations of helicity accumulation, and magnetic flux for six nonflare active 
regions. The helicity is shown as crosses, and the magnetic flux is shown as diamonds. 
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Fig. 4. — Helicity parameters with GOES X-ray flux integrated over the flaring time. Cor- 
relations of the integrated soft X-ray flux with (a) average helicity change rate of phase I, 
(6) the amount of helicity accumulation during phase I, and (c) helicity accumulation time. 
Correlation coefficient (CC) is specified in each panel. In (rf), the amount of helicity accu- 
mulation is plotted as a function of the accumulation time. The uncertainties of the average 
helicity change rate, the amount of helicity accumulation, and the helicity accumulation time 
are shown as error bars in each panel. 



